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ABSTRACT 

We discuss a possible physical connection between helium-rich (Y > 0.35) stellar populations of 
massive globular clusters (GCs) and the ultraviolet (UV) upturn of galactic spheroids by using ana- 
lytical and numerical models. In our model, all stars are initially formed as bound or unbound star 
clusters (SCs) formed from giant molecular clouds (GMCs) and the SCs can finally become GCs, open 
clusters, and field stars depending on physical properties of their host GMCs. An essential ingredient 
of the model is that helium-rich stars are formed almost purely from gas ejected from massive asymp- 
totic giant branch (AGB) stars. The helium-rich star formation is assumed to occur within massive 
SCs if the masses of the progenitor GMCs are larger than a threshold mass (-Mth res ). These massive 
SCs can finally become either massive GCs or helium-rich field stars depending on whether they are 
disintegrated or not. Using this model, we show that if the initial mass functions (IMFs) in galactic 
spheroids are mildly top-heavy, then the mass fractions of helium-rich main-sequence stars (Fn e ) can 
be as large as ~ 0.1 for Mthres = 10 7 M Q . F^ c is found to depend on IMFs and Mthr es such that it 
can be larger for shallower IMFs and smaller M t hrcs- The inner regions of galactic spheroids show 
larger Fh in almost all models. Based on these results, we suggest that if the UV upturn of elliptical 
galaxies is due to the larger fractions of helium-rich stars, then the origin can be closely associated 
with top-heavy IMFs in the galaxies. 

Subject headings: galaxies: star clusters- globular clusters: general - starsdormation - galax- 
ies: elliptical and lenticular, cD - ultraviolet: galaxies 



1. INTRODUCTION 

Recent observational and theoretical studies of the 
Galactic GCs have suggested that some of the massive 
GCs (ui Cen and NGC 2808) have significant fractions of 
helium-rich stars (e.g., Bedin et al. 2004; Norris 2004; 
Lee et al. 2005a; Piotto et al. 2005, 2007; Renzini 
2008). The observed unusual level of helium enhance- 
ment (AY/AZ w 70; Piotto et al. 2005) and sizable 
fractions of helium-rich stars in these GCs have driven 
many theoretical studies to make great efforts in under- 
standing where and how helium-rich stars can be formed 
(e.g., Bekki & Norris 2006; D'Antona & Ventura 2007; 
D'Ercole et al. 2008). As discussed by Renzini (2008), 
helium-rich stars can be formed from gaseous ejecta from 
massive AGB stars with initial masses (mi) ranging from 
3M to 8M Q , though the ejecta should not be diluted by 
interstellar gas (ISM) to keep the original high helium 
abundances of the ejecta. 

Code & Welch (1979) investigated the integrated light 
of seven elliptical and SO galaxies for a wavelength range 
of 1550-4250 A and found that some of them show an 
increase in energy at the shortest wavelengths. Since 
this discovery of the UV upturn, the origin of the UV 
upturn has been extensively discussed both in observa- 
tional and theoretical studies (e.g., Bertola et al. 1980; 
Burstein et al. 1988; Greggio & Renzini 1990; Horch 
et al. 1992; Dorman et al. 1995; Brown et al. 1997; 
O'Connell 1999; Yi et al. 1997, 1998, 2005). Although 
there can be a number of stellar candidates responsible 
for the UV upturn in galaxies, one of the most promis- 
ing ones is old horizontal-branch stars (e.g., Yi 2008 for 
a recent review). It is suggested that enhanced helium 



abundances can play an important role in the formation 
of hot stars responsible for the UV upturn (Yi 2008). 

Brown et al. (2003) investigated the UV emission in 
eight early-type galaxies at z = 0.33 (a look-back time 
of 3.9 Gyr) and found that the UV emission in these 
galaxies is significantly weaker than it is in the current 
epoch. Observational studies by Galaxy Evolution Ex- 
plorer ( GALEX) investigated the UV properties of bright 
cluster galaxies (BCGs) in clusters at z < 0.2 and com- 
pared them with those of nearby giant elliptical galaxies 
(e.g., Lee et al. 2005b; Ree et al. 2007). Reeetal. (2007) 
concluded that the observed evolution of FUV — V color 
with a model in which the dominant FUV source is hot 
horizontal-branch stars (e.g., Ree et al. 2007). Yi et al. 
(2011) investigated correlations between the strength of 
the UV upturn and global galactic properties (e.g., lu- 
minosity) in BCGs and did not find any remarkable cor- 
relations (see also Loubser & Sanchez-Blazquez 2011 for 
similar results). Yi et al. (2011) therefore concluded that 
the helium sedimentation scenario proposed by Peng & 
Nagai (2009) can not be supported by their observational 
results. 

These observational and theoretical studies appear to 
suggest that there can be helium-rich stellar populations 
in diverse objects with dramatically different masses and 
sizes ranging from GCs to BCGs. Although it would be 
likely that different astrophysical objects with helium- 
rich stars have different origins, it would be possible that 
they can have a common origin. Previous theoretical 
studies, however, have not yet provided an unified pic- 
ture for the origins of helium-rich stars in GCs, galactic 
bulges, and elliptical galaxies with the UV-upturn. It is 
therefore worthwhile to construct a model to discuss the 
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origins of helium-rich stars observed in GCs and galactic 
spheroids in a self-consistent manner. 

The purpose of this paper is to present a new model 
which can provide a possible explanation for the ori- 
gins of helium-rich stars with Y > 0.35 in GCs and 
galactic spheroids in a self-consistent manner. Based 
on the model, we mainly discuss what types of IMFs 
are required to explain the observed possible fraction of 
helium-rich stars in galactic spheroids with the UV up- 
turn (e.g., Chung et al. 2011). The most important in- 
gredient of the model is that helium-rich stars are formed 
almost purely from gas ejected from massive AGB stars 
with 3M Q <m\< 8M Q with no/little dilution of the gas 
with ISM. This formation of new stars from AGB cjccta 
have been included also in recent theoretical models of 
GC formation (e.g., D'Ercole et al. 2008). 

In the present study, we assume that such formation 
of helium-rich stars almost directly from AGB cjccta can 
be possible only within SCs formed from massive GMCs. 
By considering that the vast majority of stars are ob- 
served to be initially formed as bound or unbound SCs 
(e.g., Lada & Lada 2003), we assume that all stars form 
as SCs from GMCs in the model. Some massive SCs can 
retain gaseous ejecta from AGB stars and consequently 
have new helium-rich stars formed form the ejecta, and 
the helium-rich stars can become field stars of their host 
galaxies when the SCs become disintegrated. Such mas- 
sive SCs can finally become GCs or nuclear SC (or stellar 
galactic nuclei) depending on their birth places, if they 
are not destroyed by their host galaxies. The IMFs and 
GMCMFs are key parameter that can determine mass 
fractions of helium-rich stars in galaxies. By using this 
new model, we discuss the origin of physical properties of 
the Galactic GCs and the UV upturn in galactic bulges 
and elliptical galaxies. 

Although it remains less clear how much fractions of 
stars should be helium-rich stars in galactic spheroids 
with the UV upturn (Yi 2008), Chung et al (2011) have 
recently suggested that if about 11% of stellar popula- 
tions in galactic spheroids are helium-rich stars, then 
they can show the UV upturn for a given IMF. In their 
models, they considered that the major source of far-UV 
flux originates from metal-poor and helium-enhances hot 
horizontal branch stars. We therefore consider that if the 
fractions of helium-rich stars in galactic spheroids are less 
than ^0.1, the spheroids are unlikely to show the UV 
upturn in the present models. We investigate in what 
physical conditions (e.g., IMFs) the fractions of helium- 
rich stars can be above 0.1 for galactic spheroids to show 
clearly the UV upturn. 

A possible physical link between helium-rich stellar 
populations and the UV upturn of galactic spheroids has 
been already pointed out by a number of authors (e.g., 
Yi 2008). Thus the main point of the present study is 
not to propose the importance of helium-rich stellar pop- 
ulations in the origin of the UV upturn but to discuss 
how galactic spheroids can have significant fractions of 
helium-rich stellar populations. It should be also stressed 
that the present model is idealized and less realistic in 
some points (e.g., no chemical evolution of galaxies): the 
model should be regarded as a first step toward better 
understanding the origin of helium-rich stars in GCs and 
galaxies. More sophisticated numerical simulations in- 
cluding SC formation processes in galaxies will need to 



be done in our future studies to address the UV upturn 
problem in a much more quantitative way. 

The plan of the paper is as follows: in the next sec- 
tion, we describe a model which enables us to estimate 
(i) the mass fraction of helium-rich stars in a single GC 
and (ii) the mass fraction of helium-rich main-sequence 
(MS) stars in a galaxy for a given IMF and GMCMF. 
In §3, we present the results on the number/mass frac- 
tion of massive GCs with helium-rich stars in the Galaxy 
and the mass fractions of helium-rich MS stars in galax- 
ies with different IMFs and GMCMFs In §4, we provide 
important implications of the present results in terms of 
the origin of helium-rich stars in GCs, the Galactic bulge, 
and elliptical galaxies. We also discuss the origin of the 
radial gradients of helium-rich stars in galactic spheroids 
in this section. We summarize our conclusions in §5. 

2. THE MODEL 
2.1. A cluster disintegration scenario (CDS) 

2.1.1. SCs as fundamental building blocks of galaxy 

formation 

We adopt a scenario in which (i) all stars in galaxies are 
formed as bound or unbound SCs from GMCs and (ii) 
helium-rich stars with Y > 0.35 can be formed from gas 
ejected from AGB stars only in massive SCs. Some SCs 
can become field stars in their host galaxies after disin- 
tegration and others can become massive GCs or nuclear 
star clusters (or stellar nuclei) in galaxies. Furthermore, 
if massive SCs with helium-rich stars are destroyed by 
their host galaxies, then the stars become helium-rich 
field populations in the galaxies. Therefore, formation 
efficiencies of massive SCs and IMFs in GCs and galax- 
ies are key parameters which can determine their mass 
fractions of helium-rich stars. This scenario is referred to 
as "cluster disintegration scenario" (CDS) just for conve- 
nience in the present study. Since a number of possibly 
unfamiliar acronyms are used in the present study, their 
physical meanings are briefly summarized in Table 1 in 
order for readers to understand more clearly the present 
paper. 

2.1.2. Two stellar populations HNSs and HRSs with 

different Y 

In the CDS, stars formed directly from GMCs are as- 
sumed to have "normal" Y that is observed in HII re- 
gions of galaxies (e.g., Peimbert, Luridiana & Peimbert 
2007) and predicted from canonical chemical evolution 
(i.e., AY/AZ w 2) and they are referred to as "HNSs" 
(helium- normal stars). After these first generation of 
stars form, then the second generation of stars can form 
from gaseous ejecta of AGB stars among HNSs without 
dilution of the ejecta with ISM. These second genera- 
tion of stars can show large degrees of helium enhance- 
ment (i.e., AY/AZ > 4) which can not be expected in 
canonical chemical evolution models in which ISM and 
gaseous ejecta from stars are assumed to well mix and 
then form new stars. These second generation of stars 
are referred to as "HRSs" (helium-rich stars) just for 
convenience. As demonstrated by many recent theoret- 
ical studies (e.g., D'Antona et al. 2010), massive AGB 
stars can ejecta gas with enhanced helium abundances 
(Ay ~ 0.07). Therefore it is quite reasonable to assume 
that HRSs can have high helium abundances. Figure 1 
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illustrates a whole picture for the formation histories of 
two different populations (HNSs and HRSs) from GMCs 
in the CDS. 

Our recent numerical simulations on secondary star 
formation from AGB ejecta in SCs have shown that the 
ejecta can be converted efficiently into new stars only 
in massive SCs with masses (m sc ) larger than ~ 1O 6 M 
(Bekki 2010; 2011). This is because AGB ejecta can be 
accumulated in the deep potential wells of massive SCs 
to form high-density gaseous regions with the densities 
exceeding the threshold gas density (pthres) for star for- 
mation (pthros > 10 4 — 10 5 atoms cm -3 ). We accordingly 
consider that AGB ejecta can be converted into new stars 
without the dilution with ISM if some physical condi- 
tions are met for SCs. We assume that (i) m sc /m gmc is 
constant for all SCs and (ii) HRS formation is possible 
within SCs only if the masses (m gmc ) of their host GMCs 
can exceeds a threshold mass (Mthres)- The introduction 
of M, 



thr 



is therefore consistent with recent numerical 
simulations by D'Ercole et al. (2008) and Bekki (2011). 

2.1.3. Y in HRSs 

It depends on chemical yields of AGB stars whether 
HRSs can have Y > 0.37 expected for blue main- 
sequence (MS) stars in oj Cen and NGC 2808 (e.g., Ren- 
zini 2008) and for the UV upturn in elliptical galaxies 
(e.g., Yi 2008). Although D'Antona et al. (2005) sug- 
gested that more massive AGB stars with m s > 6M© 
can have Y = 0.40, AGB stars with lower m s (< 5M Q ) 
could have smaller Y (see Renzini 2008 for a more de- 
tailed discussion on this). It should be noted here that 
Portinari ct al. (2010) suggested helium-rich populations 
in GCs to have Y ~ 0.3 rather than Y « 0.4. Consid- 
ering these recent results, we adopt a model in which 
HRSs formed from gaseous ejecta of AGB stars with 
3M© < m s < 8M© can become helium-rich populations 
of GCs and galactic spheroids. 

As discussed in our previous study (Bekki & Norris 
2006), HRSs with higher Y can be formed from ejecta of 
Type II supernovae, if the ejecta does not mix with ISM. 
We however does not discuss whether and how HRSs 
can be formed from the ejecta of Type II supernovae in 
the present study. Some observational results suggest 
that the helium-enrichment parameter AY/AZ can be 
larger than 4 (e.g., 5.3 ± 1.4 for local K-dwarf stars of 
the Galaxy; Gennaro et al. 2010). Therefore, it would be 
possible that even in canonical chemical evolution mod- 
els, metal-rich stars with Z ~ 0.03 can have large Y 
(> 0.35) for the primordial helium content of Y = 0.24. 
Since we do not include chemical evolution of GCs and 
galaxies in the present study, we do not discuss the above 
possibility either. 

2.2. IMF 

We consider that HNSs and HRSs are formed in SCs 
with different IMFs. The adopted IMF in number is 
defined as follows: 



V'(mi) = m SCi0 mi" 



(1) 



where mi is the initial mass of each individual star and 
the slope a = 2.35 corresponds to the Salpctcr IMF. The 
normalization factor m SCi o is a function of m sc , m\ (lower 



mass cut-off), and m u (upper mass cut-off): 



x (2 - a) 



rn, 



2-a _ 



(2) 



where m\ and m u are set to be free parameters in the 
present study. The IMF slopes for HNSs and HRSs are 
denoted as ct\ and ct2 respectively, and m\ i (mi 2 ) and 
m u> i (m 1>2 ) arc for HNSs (HRSs). 

Gaseous ejecta from AGB stars with masses ranging 
from mi >ag b to m Uiag b can contribute to the formation 
of HRSs and mjAGB and m u>ag b are fixed at 3M© and 
8M Q , respectively. The total mass of AGB ejecta within 
a SC (M ag b) is accordingly described as: 



Magi 



J mi „, 



m e j'ip(m)dm, 



(3) 



™l,agb 



where m c j describes the total gas mass ejected from an 
AGB star with initial mass mi and final mass mp. We 
derive an analytic form of m c j (= mi — mp) from the 
observational data by Weidemann (2000) by using the 
least-square fitting method, and find: 



m oj = 0.916m! - 0.444. 



(4) 



Using these equations, we estimate the mass fraction of 
AGB ejecta (/ m , ag b) in each individual SC for a given set 
of IMF parameters. / m , ag b is defined as follows: 



M„ 



m , agb 



(5) 



In order to discuss an increase of helium mass in a SC 
(AMho) due to gas ejection of massive AGB stars, we 
adopt a formula used by Renzini (2008): 



AM Hc = 0.15(mi - 3)M© 



(6) 



We estimate the mass fraction of the fresh helium gas 
(/m,He) in each individual SC as follows: 



Q^AM^(m)dm 

Jm.Hc " 



(7) 



It is suggested that / m ,Hc can be 0.7 per cent for a 
Salpeter IMF and mj = O.5M (Renzini 2008). If su- 
per AGB stars with 8 < mi < 10M© can also contribute 
to the production of helium-rich gas, then / mj He can be 
0.009 (Renzini 2008). These possibly small / mj H e made 
several authors to suggest that the original GCs are much 
more massive than the present ones (e.g., Bekki & Norris 
2006). 

2.3. Mass fractions of HRSs in GCs 

We investigate the mass fractions of HRSs on the MS 
(simply referred to as "MS HRSs") in the present GCs 
(/He) In order to calculate the MS turn-off mass (mTo), 
we use the following formula (Renzini & Buzzoni 1986): 

lnm TO (£ s ) = 0.0558(logt s ) 2 - 1.3381ogi s + 7.764, (8) 

where mxo is in solar units and time t s in years. We 
assume that ages of GCs and galactic spheroids are 12 
Gyr and thus rnxo = 0.885M Q . Although the present 
results can hardly depend on age differences of <~ 300 
Myr (corresponding to the TO epoch of stars with mi = 
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3M ) between HNSs and HRSs, we consider that raxo = 
0.885M Q for HNSs and m TO = 0.889M Q for HRSs. 

The initial total mass of HRSs in a SC (m SC;2 ) is given 
as: 



m S c,2 = e s f,2-MAGB, 



(9) 



where e s f ; 2 is the star formation efficiency for stars 
formed from AGB ejecta. For a strongly bound SC to 
be formed, the SC should not lose a significant fraction 
of gas left behind from star formation: e s f t 2 needs to 
be larger than 0.5 (e.g., Hills 1980). We here adopt 
e 2 = 1.0 in the present study. The total mass of MS 
HNSs (m S c,mshns) m a SC is given as: 



^sc.mshns 



/ mip(m)dm. 

J mi,i 



(10) 



The total mass of MS HRSs (™ sc ,mshrs) in a SC as follows: 

rm-ro 

m S c,mshrs = e s f,2 / mip{m)dm, (11) 

J mi, 2 

where the normalization factor of the IMF (tp) is de- 
termined by m SCt 2- Therefore, the mass fraction of MS 
HRSs in a SC is given as: 



h 



^sc.mshrs 



He 



shrs 



(12) 



^sc,mshri 

Thus /n e depends on ai, ct2, w u ,i, mi,2, and m Ui 2- 

2.4. GMCMF 

In order to discuss the number/mass fraction of MGCs 
in the GC system (GCS) of the Galaxy and mass frac- 
tions of HRSs in galactic spheroids, we introduce GM- 
CMFs. A galactic spheroid is initially composed of nu- 
merous GMCs from which unbound and bound SCs can 
be formed. We adopt the following GMCMF for GMCs: 



*("i gmc ) = M 9i0 m gmc " 



(13) 



where /? describes the slope of a GMCMF and is observed 
to be 1.6—1.7 for galaxies in the Local Group (Rosolowski 
2005). The normalization factor M^o is a function of the 
total gas mass in a galactic spheroid (M g ), m gmC! i (lower 
mass cut-off), and m gmc , u (upper mass cut-off): 



M g x (2 - f3) 



^gmc.u ^gmc.l 



2-/3 ' 



(14) 



where m gmc j and m gmCiU are set to be 10 3 M© and 
10 8 M©, respectively, in the present study. 

We consider that ft could be different in different re- 
gions of a galaxy and in different galaxies and thus in- 
vestigate models with different (3 (1.1 < /3 < 2.0). We 
assume that star formation efficiencies (e s f.i) within host 
GMCs arc constant, though they can be different be- 
tween different GMCs. This is because we focus on mass 
fractions of HRSs in GCs and galactic spheroids and 
need to more clearly show their dependences on IMFs 
and GMCMFs. The present results do not depend so 
strongly on m gmc j and m gmc , u for reasonable ranges of 
these two. 



2.5. Formation efficiency of bound SCs in GMCs 

Some SCs get disintegrated and consequently become 
field stars in their host galaxies and others survive from 
disintegration and finally become GCs. SCs with HRSs 
can finally become massive GCs (like lo Cen and NGC 
2808), which are referred to "MGCs": other GCs with 
no HRSs are simply referred to as GCs. The total mass 
of SCs with HRSs formed in a galaxy is a key parameter 
which determines the total helium mass of the galaxy. 
Therefore, the formation efficiency of bound SCs (ebsc) 
in GMCs is one of the most important parameter in the 
present study. Here ebsc is defined such that if ebsc = 1 
(0.1), then each (one in ten) GMC can form a SC that is 
strongly bound so as to form stars from AGB ejecta. 

Although it remains observationally unclear what a 
reasonable value is for eb sc , our previous simulations 
showed that formation efficiencies of GCs in starburst- 
ing galaxy mergers can be much higher in comparison 
with isolated disk galaxies owing to rather high pressure 
of ISM (Bekki et al. 2002). The simulations suggested 
that gas can be converted into strongly bound GCs rather 
than field stars (or weekly bound SCs) in major galaxy 
mergers. These results imply that if galactic spheroids 
are formed from major mergers, they can have rather 
high 

£bsc- We discuss ebsc later in §4. 

2.6. Number and mass fractions of MGCs in the Galaxy 

We discuss number and mass fractions of MGCs 
(-Fn.mgc and F mimgc , respectively) can form "genuine 
GCs" that are observed to show O-Na anti-correlations 
and thus evidence for the presence of multiple stellar 
populations (e.g., Carretta et al. 2010) in the Galac- 
tic GCS. We consider that GMCs with m gmc > 10 6 M 
(or m sc > 1O 5 M for e = 0.1) form "genuine GCs". The 
threshold GMC mass for the genuine GCs is denoted as 
m gmc gc for convenience. We estimate F n mgc as follows: 



mgc 



N, 



(15) 



where N gc is the total number of genuine GCs and given 
as: 

e hsc ^(m)dm. (16) 

■^gmc,gc 

Likewise, -F m , mgc is estimated as follows: 



ImT C '" e s f,ie bsc m^(m)dm 



(17) 



where M gc is the total stellar mass of the genuine GCs 
and given as: 



M. 



J m 



e si .ie hsc m^{m)dm. 



(18) 



2.7. Mass fractions of HRSs in galactic spheroids 
The total stellar mass of a galactic spheroid (M s ) is 



estimated as follows: 



AL 



-L 



e s [ t im^(m)dm. 



(19) 



m gmc,l 
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The total mass of HRSs in the galactic spheroid (M Sj h rs ) 
is estimated as follows: 



s,hrs 



/>m gmCjU 

/ £ sf\ 
JM thres 



ie s f,2ebsc/m,agbm*(m)dm. (20) 



In these equations, the term ey, sc is included only in 
the equation (20), because whether SCs are bound or 
unbound does not matter in estimating the total stel- 
lar mass M s . The mass fraction of HRSs in a galactic 
spheroid (Fn e ,t) is therefore given as follows: 



^Hc,t = 



M 



s,hrs 



(21) 



Although this i<He,t enables us to understand how much 
fraction of stars can be HRSs in a galactic spheroid, it is 
different from F^ c which is more useful when the origin 
of galactic spheroids with HRSs are discussed. 

The total mass of MS stars in a galactic spheroid 
(Ms.ms) is estimated as follows: 

/•mgmc.u 

M s , ms = / e sfi i/ mS) iTO*(m)dTO, (22) 

where / ms .i is the mass fraction of MS stars among all 
stars and depends on the IMF parameters of HNSs (e.g., 
ai). The total mass of MS HRSs (M Simsrir s) is estimated 
as follows: 

M Sjmshrs = / e s f i ie s f ! 2e b sc/ms.2TO*(TO)c?TO, (23) 

JM thres 

where / ms ,2 is the mass fraction of MS HRSs among all 
stars formed as HRSs and determined by IMF parame- 
ters of HRSs. Therefore Fh is given as follows: 



M 



s,mshrs 



M, 



(24) 



In the present study, e s f,i is assumed to be indepen- 
dent of m gmc and accordingly the present results on 
i^Ho do not depend on e^i. Since we focus on Fjje 
in galactic spheroids, it is unimportant whether low- 
mass SCs are unbound (or bound) to become fields stars 
(open/globular clusters) after SC disintegration. The pa- 
rameter ebsc is therefore unimportant for GMCs with 
Ti gmc < Mthrcs- We assume that ebsc = 1 for GMCs 
that form MGCs in all models. It would be possible that 
Fu e is significantly overestimated owing to the adopted 
assumption of 6b sc = 1- We later discuss this point in §4. 
Thus free parameters are a\, a 2 , mi,i, mi, 2, m u ,i, m Ut2 , 
ft, and Mthrcs in the present study. Table 2 and 3 briefly 
summarizes the physical meanings of these parameters 
and the definition of physical quantities investigated in 
the present study (e.g., Fa e ), respectively, so that readers 
can understand more clearly the present results. 

3. RESULTS 
3.1. GCs 

Figure 2 shows how / m! H and / m ,agb depend on IMF 
parameters, a\, m^i, and m Ui i. Clearly, / mj H C is rather 
small (~ 0.005) for a canonical IMF with ot\ — 2.35, 
mi,i = 0.1M Q , and mi ;U = 100M©. This result means 
that only a small fraction of original stellar mass in a SC 
can be fresh helium gas that can be used for the forma- 
tion of HRSs for a canonical IMF. However, / m! H can 



be larger than 0.01 for top-heavy IMFs (a± < 2) if m u is 
smaller (e.g., 0.02 for ot\ and m u .i = 8M©). In extreme 
situations with m Ui i = 8M© (i.e., no Type II supernova), 
/m,He can increase as the IMF slope decreases (i.e., shal- 
lower or more top-heavy IMF). Although top-heavy IMFs 
with m Ui i = 8M Q would be highly unlikely for most GCs, 
these peculiar IMFs could be associated with the origins 
of some GCs with large fractions of HRSs. It is clear 
that / m ,He depends on such that / mj H is larger for 
larger m\ t \. 

The parameter dependences of / m ,a g b on IMF param- 
eters are essentially the same with those of /m,He- SCs 
with canonical IMFs with u\ — 2.35, m\_\ = O.lrriQ, 
and m U; i = IOOMq can have / m ,agb ~ 0.08. However 
SCs with moderately top-heavy IMFs (a ~ 2) can show 
/m,agb as large as 0.2, if = 0.5m Q and m U) i are 
smaller than 3OM (e.g., / m ,a g b ~ 0.29 for — 0.5M Q 
and for m Ui i = 8M Q ). Given that these SCs can not lose 
significant fractions of their masses owing to the IMFs 
with rather small m u .i, they are highly likely to survive 
after supernova explosions so that secondary star forma- 
tion can be possible within the SCs. If most of the AGB 
ejecta can be converted into HRSs within SCs, then the 
SCs can have significant fractions of HRSs. Thus the 
combination of moderately top-heavy IMFs and rather 
small m Uj i can significantly increase / m .agb, which sug- 
gests that SCs can contain significant fractions of HRSs 
without disintegration. The importance of such combi- 
nation of a.\ ~ 2 and low m U) i in the formation of HRSs 
was not discussed in previous studies (e.g. Bekki & Nor- 
ris 2006; Renzini 2008). 

Figure 3 shows how fn e depends on a\ for a given a 2 in 
SCs with m Ui i = O.1M , m Ui i = 100M Q , and e s f, 2 = 1. 
Clearly, SCs with canonical IMFs with a,\ = 2.35 and 
a 2 = 2.35 have a small fn e (~ 0.08). SCs with more 
top-heavy IMFs in their HNSs (i.e., smaller ai) show 
larger /h c for a given a 2 . SCs with larger a 2 (more 
"bottom-heavy") show larger /h c for a given ct\. For 
SCs with canonical IMFs in their HRSs (a 2 = 2.35) to 
have significant fractions of HRSs (/h c ~ 0.2), they need 
to have top-heavy IMFs in HNSs (ai < 2). These results 
do not depend so strongly on m\ i2 and m u . 2 , though /h c 
can be larger for smaller m u ^ 2 . If SCs have too top- 
heavy IMFs (e.g., a\ < 1.5), then they can lose most of 
their mass (more than 50%) through Type II supernova 
explosions. These SCs are likely to become disintegrated 
before secondary star formation from AGB ejecta can 
proceed. Thus it is unlikely that SCs with large /n e 
(> 0.5) can be now observed as bound GCs. 

SCs with moderately top-heavy IMFs (ai ~ 2) can 
have significant fractions of HRSs (/n e ~ 0.2), if these 
SCs can not be disintegrated due to mass loss through 
Type II supernova explosions. It should be however 
stressed that e s f.2 = 1 and no mixing of AGB ejecta 
with ISM are assumed in these investigation. Further- 
more, it is assumed that all of gaseous ejecta with AGB 
stars with different masses and lifetimes can be converted 
into new stars simultaneously. In real environments of 
GC formation, e s f,2 can be significantly smaller than 1 
and AGB ejecta could be diluted by ISM. Also, massive 
stars among HRSs can explode as supernovae and con- 
sequently prevent gaseous ejecta from AGB stars with 
lower masses from being converted into new stars. There- 
fore /ho would be likely to be lower than those estimated 
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in the present study. It should be also stressed that if 
rather small m Uj i (< 3OM ) is adopted, then /n e can be 
higher than those shown in Figure 3. 

3.2. The fraction of MGCs in the Galactic halo 

By assuming that all of the Galactic GCs were formed 
from GMCs, we here briefly discuss how much fraction of 
GCs were originally formed as MGCs with HRSs (like to 
Cen and NGC 2808). In the present study, SCs formed 
from GMCs with m gmc > Mth rC s can efficiently con- 
vert their AGB ejecta into new stars without dilution 
by ISM so that they can have HRSs. Other less mas- 
sive GC that evolve from SCs formed from GMCs with 
< m gmc < M th rcs can finally become genuine 
GCs with no/few HRSs. SCs formed from GMCs with 
m gmc < TOgmcgc can not become genuine GCs in the 
present study. Given that the Galactic GCs with signif- 
icant fractions of helium-rich populations (e.g., to Cen) 
have total masses larger than 10 6 M Q , we consider that 
-Mthrcs = 1O 7 M (for e s f,i = 0.1) is reasonable. 

Figure 4 clearly shows that although the number frac- 
tion of MGCs is small (-F^mgc ~ 0.16) for a reasonable (3 
(=1.7) and M t hros = 1O 7 M , the mass fraction (F m , mgc ) 
is quite significant (> 0.6). Both F n!mgc and F mimgc are 
larger for smaller j3 and smaller M t hres owing to the larger 
fractions of GMCs with m mgc > M trir cs- The essential 
reason for large F mmgc (> 0.5) in some models is that 
GMCMFs have slopes larger than -2 (i.e., /3 < 2). It 
should be stressed that ebsc is assumed to be constant 
for all SCs in deriving these results. If host GMCs for 
MGCs have higher ebsc than those for other GCs, then 
-Fn.mgc and F m , mgc can be even larger than those shown 
in Figure 4. 

3.3. Galactic spheroids 

Figure 5 shows that the mass fraction of HRSs (Fke,t) 
in the model with f3 = 1.7 and M t hres = 1O 7 M0 is rather 
small (~ 0.05) for canonical IMFs {a\ — 2.35). Since ebsc 
is assumed to be 1 for all GMCs with m gmc > M t hrcs in 
the present study, i<Hc,t can be significantly smaller than 
0.05 in real galaxies where ebsc can be different in differ- 
ent star-forming regions and well less than 1 in some 
star-forming regions. These results mean that galac- 
tic spheroids are unlikely to have significant fractions of 
HRSs for canonical IMFs (ai — 2.35). Clearly, F Ho ,t is 
larger for smaller f3, because there are larger fractions 
of GMCs that have m gmc > M tri rcs and thus can host 
MGCs with HRSs in galactic spheroids. Also i*H e ,t is 
larger for smaller M t hrcs for a given j3 owing to larger 
fractions of GMCs with m gmc > M tri rcs- 

The mass fractions of MS HRSs (-Fhe) can be signifi- 
cantly larger than Fa e ,t in galactic spheroids, if the mass 
fractions of MS HNSs are smaller owing to top-heavy 
IMFs. Figure 5 shows how -Fhc depends on a\ for a 
given set of 0:2, /3, and M trrr cs- As expected from Figure 
4, Fne is rather small (~ 0.05) for canonical IMFs with 
ai = 2.35 and a 2 = 2.35. Irrespective of «2, 0, and 
-Wthres, IMFs of HNSs in galactic spheroids should be 
top-heavy [a\ < 2) for the spheroids to have significant 
fractions of MS HRSs (F Ho > 0.1). It should be stressed, 
however, that the observationally suggested Fh c (~ 0.1) 
for galactic spheroids with the UV upturn (Chung et al. 
2011) can be explained by the models with mildly top- 



heavy IMFs of ai ~ 2 that are not so exotic. 

It is clear that Fh is larger for larger a 2 for a given set 
of ai, /3, and Mthres- This is because larger numbers of 
low- mass HRSs (mi < 3M ) can be formed for larger a^. 
Also, File can be larger for smaller (i and smaller M t h re s, 
because larger numbers of GMCs with m mgc > M trrr cs 
can be formed. These results mean that whether galac- 
tic spheroids can have higher Fh c (> 0.1) depends on 
physical properties of GMCs and physical processes of 
secondary star formation in MGCs. It should be stressed 
here that both e s f ;2 = 1 and eb sc = 1 are assumed in 
these estimation: top-heavy IMFs are required for galac- 
tic spheroids to have significant fractions of MS HRSs 
even in these maximum formation efficiencies of HRSs 
and MGCs. 

4. DISCUSSION 

4.1. The lost GCs with helium-rich stars 

The Galactic GCs can sink into the inner region of 
the bulge owing to dynamical friction against the halo 
within ~ 13 Gyr, if their initial masses (m gc ) is enough 
large (Binney & Tremaine 1987). By assuming that the 
Galaxy has a singular isothermal sphere, the time scale 
of dynamical friction for a GC (idf ) can be estimated as 
follows: 

<df = 2 - 3( T^ ( 2k7c } <220^<loW ^ 

where IgA, n, v c are the Coulomb logarithm, the initial 
distance of the GC from the Galactic center, and the 
circular velocity of the Galaxy. In the above estimation, 
the reference value of m gc was set to be a higher, because 
GCs can have significantly higher masses at their birth. 
The above equation means that initially massive GCs 
with HRSs have already sunk into the central region of 
the Galactic bulge owing to dynamical friction, if their n 
are less than 2 kpc. These GCs can not be observed as 
the Galactic halo GCs and contribute to the formation 
of helium-rich stars in the bulge. 

In order to estimate how much fraction of massive GCs 
in the GCS of the Galaxy has been already lost in the 
bulge, we consider that the GCS has a spherical distribu- 
tion with a density profile of p(r) oc r~ . This is consis- 
tent with that observed for the Galactic GCS (Djorgovski 
& Meylan 1994). The GCS is distributed within 35 kpc 
of the Galaxy and have a half-number radius of 5 kpc. 
We here investigate a threshold radius -Rthrcs for which 
tdf can be smaller than 13 Gyr for a given m gc . GCs with 
fi < Mthres can spiral into the center of bulge within 13 
Gyr so that they can not be observed as the Galactic halo 
GCs. It is found that ~ 49% of GCs with m gc = 1O 7 M 
have been already lost for the above adopted GCS radial 
profile (i.e., i?thrcs = 4.7 kpc). The fraction of these lost 
GCs (-Fiost) is smaller for smaller m gc : FJost is 0.37 for 
m gc = 5 x 1O 6 M and 0.06 for m gc = 1O 6 M . 

These results imply that what we can now observe as 
MGCs (e.g., lj Cen and NGC 2808) could be [51-63]% of 
original MGCs (m gc > 5 x 1O 6 M ) formed in the Galac- 
tic halo. The lost MGCs would have been destroyed by 
the Galactic bulge and consequently their HRSs would 
have been dispersed into the inner bulge region: the stars 
could be now observed as HRSs in the bulge. This se- 
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lective loss of MGCs due to dynamical friction might 
well occur in galactic bulges and giant elliptical galaxies. 
Sohn et al. (2006) revealed a large number of GCs with 
strong UV light in M87 (see also Kaviraj et al. 2007). 
If these GCs contain large fractions of HRSs, then M87 
could have already swallowed a fraction of these possible 
MGCs to add their HRSs to its main stellar spheroid. 

4.2. Top-heavy IMFs for the origin of the UV upturn in 
elliptical galaxies and BCGs 

The present study has demonstrated that if IMFs are 
top-heavy (i.e., a.\ < 2), then significant fractions (~ 0.1) 
of MS stars in galactic spheroids can be HRSs. This is 
mainly because a larger amount of AGB ejecta with high 
helium abundances can be converted into new stars with- 
out dilution of the ejecta with ISM. This therefore implies 
that if the UV upturn is due largely to large fractions 
of MS HRSs within them, as observationally suggested 
(Chung et al. 2011), then galactic spheroids with mildly 
top-heavy IMFs are more likely to show the UV upturn. 
Then are the proposed top-heavy IMFs consistent with 
other observational properties of galactic spheroids such 
as their luminosity evolution and chemical abundances ? 

A recent observation suggested that IMFs in massive 
galaxies at 0.02 < z < 0.83 are significantly flatter (cor- 
responding to a = 1.3 in our model thus even shallower 
than the required slope of ~ —2) than the present-day 
value of the Galaxy and proposed a "bottom-light" IMF 
for massive galaxies at higher redshifts (van Dokkum 
2008). This observational result is consistent with the 
present proposal and accordingly implies that if the ob- 
served massive galaxies in van Dokkum (2008) are pro- 
genitors for the present galactic spheroids, then some 
of them can show the UV upturn owing to larger Fh c - 
It should be stressed here that the observed fraction of 
spheroid populations with the UV upturn is small (e.g., 
Yi et al. 2011) and thus the observed mean properties of 
high redshift spheroids could be different from those of 
the present ones with the UV-upturn. 

The required top-heavy IMFs are also consistent with 
the IMFs proposed for explaining the chemical proper- 
ties of galactic spheroids (e.g., bulges), such as metallic- 
ity distribution functions (MDFs) and [a/Fe] (e.g., Mat- 
teucci & Brocato 1990; Nagashima et al. 2005; Ballero 
et al. 2007; Tsujimoto et al. 2010). Loewenstein (2006) 
showed that rapid star formation with top-heavy IMFs 
in massive galaxies are necessary to explain the observed 
Fe abundances of intra-cluster medium (ICM) which con- 
tains heavy metals from supernovae of massive early- 
type galaxies in cluster of galaxies. The top-heavy IMFs 
adopted in these previous studies suggest that the re- 
quired IMF for explaining the UV-upturn observed in 
some spheroids is neither unreasonable nor unrealistic. 
Larson (1998) suggested that gas clouds with higher tem- 
perature in galaxies at high redshifts can be responsible 
for star formation with top-heavy IMFs. Thus the re- 
quired top-heavy IMFs could be closely associated with 
high-redshift formation of the vast majority of stars in 
galactic spheroids. 

It should be also stressed that if IMFs are too top- 
heavy (e.g., ct\ < 1.5), then the SCs could become dis- 
integrated before AGB ejecta can be recycled and con- 
verted into new stars. About 45% in mass can become 
stars with masses larger than 8M (thus supernova) in 



SCs for a-i = 1.95, m u = 0.1M o , and m Uj i = 100M Q . 
These SCs can not become disintegrated and thus can 
continue secondary star formation, because more than 
50% of their masses can still remain in SCs (e.g., Hills 
1980). However, if ct\ = 1.5, m\ t \ = 0.1M Q , and 
m u .i = 100M Q then the mass fractions of supernovae in 
SCs are 0.74 so that SCs can lose most of their original 
masses. These SCs are highly likely to get disintegrated 
shortly after supernova explosions owing to their mass 
loss. Therefore, too top-heavy IMFs (ai < 1.5) are not 
ideal for galactic spheroids to have helium-rich stellar 
populations. 

Peng & Nagai (2009) proposed that sedimentation 
of helium in clusters of galaxies can be responsible 
for the formation of HRSs in BCGs. They also pre- 
dicted that the UV flux strength is stronger in more 
massive, low-redshift, and dynamically relaxed BCGs. 
Their predictions, however, have been recently chal- 
lenged by an observational study by Yi et al. (2011) 
which have found no correlation between the UV strength 
and rank/luminosities of BCGs and showed no clear dif- 
ference in UV upturn fraction or strength in BCGs. The 
present study provides an alternative explanation for the 
origin of the BCGs with the UV upturn: BCGs show 
the UV upturn because IMFs at their formation were 
top-heavy (ai < 2) and therefore they have larger Fn e . 
In the CDS with top-heavy IMFs, some BCGs can not 
show the UV upturn, because most stars are formed as 
SCs with canonical IMFs (whereas some can owing to 
top-heavy IMFs). 

However, the origin of UV upturn in BCGs would not 
be so simple as the CDS with top-heavy IMFs envis- 
ages. Recent theoretical models have shown that hierar- 
chical merging of smaller galaxies can play a vital role in 
the formation of BCGs (e.g., De Lucia & Blaizot 2007). 
The IMFs are observationally suggested to be different 
between faint and luminous galaxies (e.g., Hoversten & 
Glazebrook 2008). Therefore, if BCGs were formed by 
numerous mergers between galaxies with different lumi- 
nosities in clusters, then not just IMFs but also merging 
histories (e.g., fractions of faint/luminous galaxies) can 
be key determinants for whether BCGs can show the UV 
upturn. 

4.3. Halo-spheroid connection 

Recently Martell et al. (2011) have revealed that about 
3% of the Galactic halo stars can have depleted carbon 
and enhanced nitrogen abundances that are very similar 
to the chemical abundances observed for the "second- 
generation" of stars in GCs. Since these stars with 
characteristic chemical abundances can be formed from 
gaseous AGB ejecta of the "first generation of stars" 
(e.g., Bekki et al. 2007; D'Ercole et al. 2008), they sug- 
gested that (i) these stars originate from GCs and (ii) a 
minimum 17% of the present-day mass of the Galactic 
stellar halo was originally formed in GCs. Although it is 
difficult for observational studies to directly estimate he- 
lium abundances of stars in GCs, Bragaglia et al. (2010) 
have recently estimated the average enhancement in the 
helium mass fraction Y between the first and second gen- 
eration (corresponding to HNSs and HRSs, respectively, 
in the present study) for 19 GCs with the Na-0 anticor- 
relation. The estimated average enhancement is about 
0.05 — 0.11, which means that a significant fraction of 
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star in these GCs can have HRSs. If - 10% of all GC 
stars are HRSs, then about 1.7% of the halo stars can be 
regarded as HRSs. 

This smaller number of 1.7% implies that the Galac- 
tic halo can not be identified as a galactic component 
with the UV upturn. Then what mechanism is respon- 
sible for the possible large difference in the mass frac- 
tion of HRSs (i*Hc) between the Galactic stellar halo 
and galactic spheroids ? The present CDS suggests that 
most stars in the Galactic halo can originate from disin- 
tegration of low-mass SCs with canonical IMFs whereas 
those of galactic spheroids originate from stars initially in 
more massive SCs with top-heavy IMFs. Furthermore, as 
discussed in §4.1, the more massive SCs (or GCs) with 
HRSs can rapidly sink into the central regions of their 
hosts so that they can not contribute to the formation 
of HRSs in their halo regions: these more massive GCs 
can preferentially become the building blocks of the cen- 
tral spheroidal components in galaxies. Thus it is highly 
likely that Fu e can be significantly different between ha- 
los and spheroids in galaxies. 

4.4. High formation efficiencies of bound massive SCs 
at the epoch of spheroidal formation 

It should be stressed here that eb sc ~ 1 has been so 
far assumed: GMCs with m gmc > M t hros need to almost 
always form bound SCs that can finally form helium-rich 
stars. If e b sc ~ 0.1, then it is very hard for galactic 
spheroids to have Fh c ~ 0.1 even for very top-heavy 
IMFs (e.g., a\ = 1). Only spheroids with a\ ~ 1 and 
m Uj i ~ 8M Q (i.e., peculiar top-heavy IMFs) can have 
Fa c as large as 0.1 for e i, sc ~ 0.1. Photometric studies 
of super star clusters (SSCs) in a starbursting luminous 
infrared galaxy Arp 220 showed that the nuclear SSCs 
in Arp 220 contribute to ~ 20% of the total bolomct- 
ric luminosity of Arc 220 (Shioya et al. 2001). Larsen 
& Richtler (2000) found that the formation efficiency of 
SCs can be higher in local regions with high star forma- 
tion rate per unit area for disk galaxies. These obser- 
vations imply that (i) a significant fraction of new stars 
are formed as SCs in a starburst or in high-density star- 
forming regions and thus (ii) galactic spheroids can have 
high cluster formation efficiencies if they were formed 
from massive starbursts in high-density environments. 

However, owing to the lack of extensive observational 
studies on formation efficiencies of massive SCs in galax- 
ies, it is currently very hard to discuss what a reasonable 
value of e t, sc is for massive GMCs that can be progen- 
itors for bound SCs with HRSs. If e b sc can be signifi- 
cantly larger than 0.1, then the present CDS is promis- 
ing as the origin of the UV upturn in galactic spheroids, 
though top-heavy IMFs are still required. Thus, ulti- 
mately speaking, the origin of the UV upturn is closely 
related to the formation processes of MGCs from massive 
GMCs at the epoch of galaxy formation in the present 
CDS. 

4.5. Helium-rich populations in the Galactic bulge ? 

Nataf et al. (2011) revealed that the metal-rich stel- 
lar populations of the Galactic bulge can have Y ~ 0.35 
based on the observational results on physical properties 
of the red giant branch bump (RGBB). They provided 
some important implications of the metal-rich HRSs in 



the bulge and suggested that galactic bulges in general 
can have HRSs like the bulge. The present results imply 
that if the Galactic bulge had a top-heavy IMF at its 
formation, then the present stellar populations can have 
a significant fraction (~ 0.1) of HRSs. Given that previ- 
ous chemical evolution models of the Galactic bulge (e.g., 
Ballcro et al. 2007; Tsujimoto et al. 2010) suggested a 
moderately top-heavy IMF (a ~ 2.05) for the bulge, it 
is possible that the Galactic bulge can have a significant 
fraction of HRSs owing to the top-heavy IMF. 

As suggested by Nataf et al. (2011), the dominant 
populations of the RGBB (with [Fe/H]> 0) in the bulge 
can have a rather high Y (~ 0.35): not just an ~ 10% 
of the bulge population needs to have Y ~ 0.35 to ex- 
plain their observational results. If more than 30% of 
the metal-rich stellar populations have Y ~ 0.35, then 
the IMF slope (a) for the populations should be well 
less than 1.5 for (3 = 1.7 and Af t hrcs = 10 7 M Q . The 
required IMF is significantly more top-heavy than those 
suggested by chemical evolution studies (e.g., a ~ 2.05; 
Tsujimoto et al. 2010) and thus would not explain the 
observed global properties of the bulge such as the metal- 
licity distribution function of the bulge stars. However, 
if metal-rich stellar populations of the bulge has a very 
top-heavy IMF and if other populations have moderately 
top-heavy IMFs, then it would be possible that both 
the observed high fraction of metal-rich HRSs and global 
chemical properties in the bulge can be self-consistently 
explained. This possibility needs to be explored in our 
future studies based on detailed chemical evolution mod- 
els. 

4.6. The origin of the correlation between the strength 
of UV upturn and Mg2 index 

Recent observational studies of 48 nearby early-type 
galaxies by the SAURON project (Bureau et al. 2011) 
have confirmed the presence of a negative correlation be- 
tween FUV - V color and Mg line strength originally 
proposed by Burstein et al. (1988). The present study 
provides the following possible explanation for the origin 
of this "Burstein relation" (i.e., a correlation of Mg2 with 
(1550 — V) color). More massive elliptical galaxies can 
retain more efficiently the ejecta of supernovae so that 
they can finally have higher metallicities of their stars 
(e.g., Arimoto & Yoshii 1987). If more massive elliptical 
galaxies have shallower (i.e., more top-heavy) IMFs, then 
they can have large Fu e and thus show the stronger UV 
upturn. Therefore more massive (or luminous) elliptical 
galaxies can have higher metallicities and thus higher 
Mg2 as well as the stronger UV upturn. Given that 
the observed positive correlations between luminosities, 
Mg indices, and velocity dispersions in elliptical galaxies 
(e.g., Faber & Jackon 1976; Bender et al. 1992), more 
massive (or luminous) elliptical galaxies can have higher 
Mg2, higher stellar velocity dispersion, and the stronger 
UV upturn. Thus, the dependence of IMF slopes on 
galactic masses/luminosities can be responsible for the 
origin of the Burstein relation. 

Although Hoversten & Glazebrook (2008) revealed 
that galaxies significantly fainter than the Galaxy show 
steeper IMFs, it remain observationally unclear how IMF 
slopes depend on global galactic properties and forma- 
tion environments of galaxies. Therefore it is too prema- 
ture to conclude whether the CDS with IMF variation is 
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promising or not. As reviewed by Yi (2008), only a small 
fraction of elliptical galaxies show the UV upturn, which 
needs an explanation. If elliptical galaxies are formed 
from merging of smaller galaxies with different IMFs, 
then the origin of the UV upturn would not be so simple 
as the CDS explains. 

The observed strong correlation between Mg2 and 
(1550 — V) color (Burstein et al. 1988) suggests that 
metallicities play a role in the formation of the correla- 
tion. Figure 1 in Burstein et al. (1988) also showed a 
large difference in (1550 — V) colors (~ 2.5 mag) among 
elliptical with different Mg2 ranging from <~ 0.20 to 0.36. 
It is unclear whether this large difference can be quan- 
titatively explained by IMF variation alone, because the 
present study can not predict (1550 — V) colors as a func- 
tion of Fh c - It would be possible that the combination 
of high metallicities and large Fh can make (1550 — V) 
colors significantly bluer in elliptical galaxies. It is ac- 
cordingly our future study to include chemical evolution 
and thereby to discuss the dependences of Fn e on metal- 
licities in elliptical galaxies. 

4.7. Radial gradients of helium-rich stars in galactic 
spheroids 

Carter et al. (2011) have recently investigated radial 
gradients of the FUV excess in 52 galaxies observed by 
GALEX and found that some of them show a positive 
gradient in the (FUV - NUV) color. They therefore sug- 
gested that the observed gradient can be due to a helium 
abundance gradient and furthermore that the presence 
of the gradient can place a strong limit on the impor- 
tance of dry mergers in elliptical galaxy formation. In 
order to discuss the origin of the observed radial gra- 
dients of (FUV - NUV) color, we have constructed a 
toy N-body model in which a galactic spheroid is formed 
from merging of numerous SCs with the mass fraction of 
HRSs being a free parameter. The details of the model 
are given in the Appendix A. 

Figure 7 shows time evolution of spatial distributions 
of HNSs and HRSs in a galactic spheroids composed ini- 
tially of numerous SCs for the standard model. As SCs 
closely interact with one another in early evolution phase 
of the galactic spheroid, smaller and less massive SCs are 
destroyed by the tidal field of the host spheroid and by 
larger and more massive SCs. The stars (HNSs) initially 
in less massive SCs are dispersed during destruction of 
their host SCs and finally become field HNSs. On the 
other hand, massive SCs with m gmc > Mthres and thus 
with HRSs can not be destroyed efficiently during early 
multiple SC interaction. These massive SCs can sink into 
the inner region of the spheroid owing to dynamical fric- 
tion against field stars (HNSs) and then interact dynam- 
ically will other massive SCs there. Consequently, the 
massive SCs can get disintegrated there and their HRSs 
can be dispersed to become field HRSs there. Since the 
inner region of the spheroid can be finally dominated by 
stars originally from HRSs of massive SCs, the inner re- 
gion can have a larger fraction of HRSs within ~ 1 Gyr 
dynamical evolution. 

Figure 8 shows how radial -FHe,t gradients depend on 
models parameters, M t hres and 7, for a give f3 (=1.7). 
Clearly, all models shows negative gradients in the sense 
that inner regions of galactic spheroids have higher frac- 
tions of HRSs. For example, -Fko^t can be ~ 0.08 within 



the central 200 pc in the standard model, though it is 
rather small (~ 0.01) in the outer region (R ~ 1.8 kpc). 
The parameter Mthres can control the central value of 
i<He,t such that -Fne.t can be larger for smaller M t hres- 
The i^He.t gradients can be smaller in models with 7 = 
in which there is no mass-size relation for SCs. The de- 
rived rather flat radial distribution is due to higher de- 
grees of dynamical mixing of SCs with different mass in 
these models. Since Fh c is proportional to Fn ,t for a 
given set of IMFs, these results on Fn e ,t gradients can 
be true for i*He- The present study accordingly demon- 
strates that galactic spheroids can show negative radial 
gradients of i^He.t and Fn e and therefore negative gradi- 
ents of the strength of the UV upturn (i.e., stronger in 
their inner regions). The present study thus can clearly 
explain the above observational result by Carter et al. 
(2011), if the observed color gradients are due largely to 
helium abundance gradients. 

5. CONCLUSIONS 

We have adopted a model in which (i) all stars are 
formed as bound or unbound SCs from GMCs and (ii) 
HRSs with Y > 0.35 can be formed from gas ejected from 
AGB stars only in massive SCs with m gmc > M t hres- 
Some SCs can become field stars in their host galaxies 
after disintegration and others can become massive GCs 
or nuclear star clusters (or stellar nuclei) in galaxies. 
If massive SCs with HRSs are destroyed by their host 
galaxies, then HRSs become helium-rich field stars in 
the galaxies. We have investigated mass fractions of 
MS HRSs (/n e and Fh c ) in GCs and galactic spheroids 
and their dependences on IMFs, M t hrcs, and GMCMFs 
based on the CDS ("cluster disintegration scenario"). 
We have also investigated the radial gradients of Fn e 
in galaxies with mass-size relations of GMCs by using 
N-body simulations on dynamical evolution of multiple 
cluster systems. We summarize our principle results as 
follows. 

(1) The mass ratios (/ m ,agb) of gaseous ejecta from 
massive AGB stars with 3M Q < m\ < 8M Q to initial 
total masses of SCs are typically ~ 0.08 for a canonical 
IMF (cti = 2.35) with the reasonable lower and upper 
cut-off masses. If these gas can be efficiently converted 
into new stars with enhanced helium abundances (i.e., 
HRSs) in SCs, then the mass fractions of MS HRSs 
among all MS stars (/h c ) are <~ 0.08 for canonical IMFs 
and a SC age of 12 Gyr. Only if original SCs have 
top-heavy IMFs (e.g., ct\ ~ 1.5), then /n e can be larger 
than 0.3 'for a 2 > 1.85. 

(2) The initial number and mass fractions of massive 
GCs (MGCs) with HRSs (-F n ,mgc and .F mingc , respec- 
tively) among the Galactic GCs formed from GMCs 
with ra gmc > 10 6 M Q is ~ 0.16 and 0.66, respectively, 
for M th = 1O 7 M and /3 = 1.7. Both F n , mgc and F miingc 
are lager for smaller M tri and smaller (3. MGCs with 
m gc = 5 x 1O 6 M can sink into the Galactic center owing 
to dynamical friction to disappear from the halo region 
within 13 Gyr, if the initial radii from the Galactic 
center are less than 3.4 kpc. About 37% (49%) of the 
Galactic GCs with m gc > 5 x 1O 6 M (m gc > 1O 7 M ) 
have already been lost from the halo owing to efficient 
dynamical friction of these massive GCs. Thus, the 
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Galactic halo initially could have a larger fraction of 
MGCs. 

(3) The mass fractions of HRSs among all stars (-Fke.t) 
in galactic spheroids are ~ 0.05 for M t hres = 10 7 M Q , 
fi = 1.7, / m ,agb = 0.1, and e bsc = 1. i*k e ,t of galactic 
spheroids are larger for smaller M t h and smaller j3 for 
a given IMF. The mass fractions of MS HRSs among 
all MS stars (-Fke) m galactic spheroids are 0.06 for 
canonical IMFs \a x = a 2 = 2.35), M th = 1O 7 M , 
P = 1.7, and a galaxy age of ~ 12 Gyr. Ffjc in galactic 
spheroids can be larger for smaller M t h and smaller /? 
for a given IMF, and it can be significant (> 0.1), if 
original SCs (i.e., building blocks of the spheroids) have 
top-heavy IMF (e.g. a\ < 2). These results suggest that 
if the observed UV upturn in bright elliptical galaxies 
is due to the larger fractions (<~ 10%) of HRSs within 
them, the IMFs need to be top-heavy. If eb sc is as small 
as ~ 0.1, then bright elliptical galaxies are unlikely to 
show the UV upturn owing to the small fractions of 
HRSs. 

(4) The Galactic bulge can have a larger Fn c if it had 
a top-heavy IMF at its formation in the CDS. Given 
a number of suggestions on the top-heavy IMF of the 
bulge by chemical evolution models of the bulge, the 
observed possible presence of helium-rich stars with 
Y ~ 0.35 in the bulge (e.g., Nataf et al. 2011) can be 
understood in the context of the CDS. HRSs in MGCs 
that were initially located in the Galactic halo and had 
sunk into the bulge can be currently observed as field 
HRSs in the bulge, though the contribution of such stars 
to the entire helium-rich population is rather minor. It 
is doubtlessly worthwhile for observational studies to 
investigate whether or not the bulge has a larger fraction 
of HRSs in the inner regions. 

(5) If Mth and (3 are constant in elliptical galaxies with 
different physical properties, then elliptical galaxies with 
shallower (or more top-heavy) IMFs can show larger 
fraction of HRSs. Therefore, if more massive/luminous 
elliptical galaxies have shallower IMFs, then they can 
have larger fraction of HRSs. More massive elliptical 
galaxies can retain more efficiently the ejecta of super- 
novae so that they can finally have higher metallicities 
of their stars (e.g., Arimoto & Yoshii 1987). Therefore, 
if the Burstein relation (a correlation between UV flux 
and Mg index; Burstein et al. 1988) in bright elliptical 



galaxies is due to the correlation between Fn c and Mg 
indices, then the relation can be understood in terms 
of shallower IMF slopes (ai) in more massive/luminous 
elliptical galaxies. 

(6) Galactic spheroids can show negative radial 
gradients of -Fke.t & n d Fn e (i.e., higher in inner parts) 
for a reasonable set of model parameters. This is 
mainly because massive SCs with HRSs can sink rapidly 
into the inner regions of galaxies owing to dynamical 
friction and disintegrate there. The HRSs initially in 
the SCs can be dispersed to become field HRSs there 
after disintegration of the host SCs. Therefore they can 
show the stronger UV upturn in their inner regions, 
if Fue can determine the strengths of the UV upturn. 
The observed spatial distribution of the FUV excess 
in elliptical galaxies can be understood in the context 
of formation and evolution of MGCs with HRSs. The 
central regions of galactic spheroids can be composed 
of two different stellar populations with normal helium 
abundances and enhanced ones (i.e., HNSs and HRSs). 

Thus, an advantage of the CDS is that both (i) the ab- 
sence or presence of the UV upturn in galactic spheroids 
and (ii) correlations between the strength of the UV up- 
turn and galactic properties (e.g., Mg2 index) can be 
discussed in the context of IMF slopes of galaxies in a 
self-consistent manner. A disadvantage of the CDS is 
that it is unclear whether and in what physical condi- 
tions ebsc can be ~ 1 for massive GMCs in galaxies. The 
formation of HRSs is possible, if helium-rich gas from 
AGB stars is converted into new stars without dilution 
of the gas by ISM with normal helium abundances. In the 
CDS, the inner regions of massive SCs are assumed to be 
the production sites of HRSs. It would be possible that 
AGB ejecta can be converted into new cold gas clouds 
outside SCs and consequently into new stars in deep po- 
tential wells of galactic spheroids, if there is no/little cold 
gas with normal helium abundances there. We need to 
investigate where and how AGB ejecta from field stars 
and SCs can be converted into new stars in galaxies us- 
ing a more sophisticated chemodynamical simulation in 
our future studies. 
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APPENDIX 

N-BODY MODELS FOR MULTIPLE SC EVOLUTION 

We investigate dynamical evolution of a galactic spheroid composed of numerous SCs using collisionless N-body 
simulations based on an idealized model for galactic spheroids. The main aim of this numerical investigation is to 
illustrate the radial gradients of the mass fractions of HRSs (Fn c t (R)) in the spheroids. This idealized model is used 
only in the present preliminary study for the origin of helium-rich stars in galactic spheroids, and more sophisticated 
models including chemical evolution due to gaseous ejection from supernovae and AGB stars will be used in our future 
studies. However, we consider that this "toy" model enables us to grasp some essential ingredients of the formation 
and evolution of radial -Fke.t gradients in galactic spheroids. 

A galactic spheroid has an initial total mass M ga i and a size i? ga i and composed of SCs with the total number of 
iV SCit . The radial distribution of SCs are described by a Plummber profile (e.g., Binney & Tremaine 1987) with the 
scale-length of 0.2i? ga i. Each SC with a mass m sc and a size r sc is composed of many stars and has a Plummber density 
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profile with the scale-length of 0.2r sc . A galactic spheroid is assumed to be initially in dynamical equilibrium, so the 
velocities of each SCs is given according to the velocity dispersion profile of the adopted Plummer model. 

A SC is composed of HNSs and HRSs and these two populations have different initial distributions within the SC. 
Recent numerical simulations have shown that secondary star formation from AGB ejecta can proceed mostly in the 
central regions of MGCs (e.g., D'Ercole et al. 2008; Bekki 2010, 2011). Guided by these results, HNSs and HRSs are 
located at 0.1r sc < R < r sc and R < 0.1r sc , respectively. The mass fraction of HRSs are calculated for a given /3 and 
Mthres- The total particle number in a SC is proportional to m sc so that masses of stellar particles can be all the same. 

A progenitor GMC for a bound/unbound SC has a mass m gmc and a size r gmc . If we use the observed relation 
between mass densities and sizes of GMCs discovered by Larson (1981) and the observed typical mass and size of 
GMCs in the Galaxy (e.g., Solomon ct al. 1979), then m gmc — r gmc relation can be described as follows: 

where 7 ~ 0.5. We investigate models with different 7, because initial SCs can have mass-size relations different from 
those of GMCs. We consider that when SCs are formed, SCs have a mass-size relation similar to the above one and 
thus assume that r sc can be the same as r gmc described above for a given m gmc . 

We mainly present the results of the "standard" dynamical models in which M ga \ = 1O 9 M0, i? ga i = 2 kpc, /? = 1.7, 
-Mthres = 1O 7 M m gmCi i = 10 4 M Q , m gmCiU = 10 8 M Q , and 7 = 0.5, though we investigate models with different 
parameters. In order to perform numerical simulations, we use the latest version of GRAPE (GRavity PipE, GRAPE- 
DR), which is the special-purpose computer for gravitational dynamics (Sugimoto et al. 1990). The total particle 
number used in a simulation is ~ 200000, which we consider to be reasonable for this preliminary investigation. The 
gravitational softening length is set to be 20 pc for all models. We focus only on the final radial gradients of i<Hc,t in 
galactic spheroids formed from dynamical evolution of numerous SCs. 
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TABLE 1 

Meaning of acronyms 



Acronyms Meanings 



CDS Cluster disintegration scenario 

SC Star cluster 

GC Globular cluster 

HNS Helium-normal stars 

HRS Helium-rich stars 

MS Main-sequence (stars) 

MGC Massive GCs (i.e., SCs with HRSs) 

GMCMF CMC mass function 



TABLE 2 

Model parameters and their physical meanings 



Parameters Physical meanings 



ai The IMF slope for HNSs 

a 2 The IMF slope for HRSs 

mi , 1 The lower mass cut-off of the IMF for HNSs 

m u , 1 The upper mass cut-off of the IMF for HNSs 

mi , 2 The lower mass cut-off of the IMF for HRSs 

m u , 2 The upper mass cut-off of the IMF for HRSs 

/3 The slope of the GMC mass function (GMCMF) 

M t h r es The threshold mass for MGC formation 

ebsc The formation efficiency of bound SCs within GMCs (fixed) 



TABLE 3 

The definition of physical quantities investigated in the present study 



Quantity Meaning 



/ m ,Hc The mass fraction of fresh helium gas in a SC 

/m.agb The mass fraction of AGB cjecta in a SC 

jf He The mass fraction of MS HRSs in a SC 

Fn,mgc The number fraction of MGCs with HRSs in all genuine GCs 

F m ,mgc The mass fraction of MGCs with HRSs in all genuine GCs 

F ti H e The mass fraction of HRSs in a spheroid 

i<He The mass fraction of MS HRSs in a spheroid 
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He-normal/rich SCs 



Fig. 1. — An illustration for the whole picture of the star formation processes from GMCs in the CDS. The evolution of newly formed 
HNSs (blue stars) depends on the masses of their host GMCs (m gmc ). HNSs formed in lower mass GMCs with TTigmc ^ ^^thres can not form 
HRSs (red stars) from gaseous ejecta of their AGB stars owing to the shallow gravitational potentials of the SCs. The HNSs can finally 
evolve either field stars after disintegration of their host SCs or into bound SCs with only HNSs. On the other hands, if m gmc > Merest 
then gaseous ejecta of AGB stars among HNSs can be converted into new HRSs (red stars) owing to deeper gravitational potentials of 
the SCs. The stars in the SCs can become either field HNSs/HRSs or main components of massive SCs with HRSs depending on physical 
properties of the SCs (i.e., on whether the SCs become disintegrated). The field HRSs in galactic spheroids are responsible for the UV 
upturn in the CDS. 
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Fig. 2. — The dependences of the mass fractions of fresh helium gas (/ m ,He; upper) and those of the AGB ejecta (/ mjag b; lower) 
in SCs on the IMF slopes of SCs (cti) for different lower (mj i) and upper (m Uj i) mass cut-offs: mi i = O.IMq and m u ,i = IOOMq 
(blue), m\ i = 0.5Mq and m Ui i = IOOMq (red), mn = O.IMq and m Uj i = 30Mq (cyan), m\ i = 0.5Mq and m Uj i = 30Mq (magenta), 
mi ! = O.IMq and m Uj i = 8Mq (green), and m\ i = 0.5Mq and m U; i = 8Mq (black). Here ot\ = 2.35 corresponds to the Salpctcr IMF. 
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Fig. 3. — The dependences of mass fractions of MS HRSs (/h c ) on a i m SCs for mj 2 = O.IMq and m Uj 2 = 8Mq (upper), and 
m h2 = O.IMq and m u , 2 = 1OOM (lower) for different IMF slopes for HRSs (a 2 ): a 2 = 1.35 (blue), a 2 = 1.85 (red), a 2 = 2.35 (cyan), 
02 = 2.85 (magenta), and a 2 = 3.35 (green). Here a\ = 2.35 corresponds to the Salpeter IMF. 
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Fig. 4. — The dependences of number and mass fractions of massive GCs with HRSs (-F n ,mgc and -F m!mgc , respectively) in the Galactic 
"genuine" GCs (formed from GMCs with m gmc > 1O 6 M ) on GMCMF slopes (/3) for different threshold masses for MGC formation 
(M thres ): M thres = 3 X 1O 6 M (blue), M t hv CS = 1 X 10 7 M Q (red), Af thrcs = 3 X 10 7 M Q (cyan), and Mthrcs = 5 X 1O 7 M (magenta). 




1.2 1.4 1.6 1.8 

P 

Fig. 5. — The dependences of the mass fractions of HRSs (Jjiet) m galactic spheroids on GMCMF slopes (/?) for canonical IMFs 
{at = 2.35) and different threshold masses for MGC formation (M thros ): M thros = 3 X 10 6 M Q (blue), M thrcs = 1 X f0 7 M Q (red), 
A^thres = 3 X 1O 7 M (cyan), and M thrcs = 5 X 1O 7 M (magenta). 



18 




Fig. 6. — The dependences of mass fractions of MS HRSs in galactic spheroids (-Fke) on ot\ in the models with = 1.4 and M t ^ ICS = 1O 7 M0 
(top), /3 = 1.7 and M t i lros = 10 7 Mq (the second from top), /3 = 1.4 and M t h ros = 3 X 1O 7 M0 (the second from bottom), /3 = 1.7 and 
-^thres = 3 X 1O 7 M0 (bottom) for different ol^'- &2 = 1-35 (blue), «2 = 1-85 (red), cti = 2.35 (cyan), ot^ = 2.85 (magenta), and cti = 3.35 
(green) . 
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Fig. 7. — Time evolution of the distributions of HNSs (blue) and HRSs (red) projected onto the x-y plane for the standard model. The 
time (T in units of Gyr) that has elapsed since the simulation starts is shown in the upper left corner of each panel. 
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Fig. 8. — The radial gradients of mass fractions of HRSs (i<He,t) in models with different parameters: 7 = 0.5 and M t hrcs = 1O 7 M0 
(blue), 7 = 0.5 and M thlcB = 3 X 10 7 Mq (red), 7 = 0.5 and M thrcs = 5 X 10 7 Mq (green), and 7 = 0.0 and M thrcB = 10 7 Mq (magenta). 



